Novel Electrostatic Attraction from Plasmon Fluctuations 
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In this Letter, we show that at low temperatures, zero-point fluctuations of the plasmon modes 
of two mutually coupled 2-D planar Wigner crystals give rise to a novel long-range attractive force. 
For the case where the distance d between two planar surfaces is large, this attractive force has an 
unusual power-law decay, which scales as d" 7 ^ 2 , unlike other fluctuation-induced forces. Specifically, 
we note that its range is longer than the "standard" zero-temperature van der Waals interaction. 
This result may in principle be observed in bilayer electronic systems and provides insight into the 
nature of correlation effects for highly charged surfaces. 
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Fluctuation-induced forces are ubiquitous in nature 
(for a recent review, see Q) and constitute an important 
contribution to the interactions of many condensed mat- 
ter systems . The classic example is the Casimir effect 
H in which quantum fluctuations of the electromagnetic 
field between two parallel conducting plates lead to an 
attractive force between them. In the context of statis- 
tical physics, Fisher and de Gennes Q have suggested 
that a similar effect also exists at or near the critical 
point of a system which is confined between two planes. 
Other examples include colloid particles immersed in a 
critical fluid ||, superfluid films (6j, liquid crystals |Q, 
and protein inclusions in fluctuating membranes In 
general, fluctuation-induced forces arise because exter- 
nal constraints modify the fluctuations of a correlated 
medium. These interactions, which are usually long- 
ranged, are controlled by thermal fluctuations at finite 
temperature or quantum fluctuations at low tempera- 
ture. In this Letter, we present arguments for a long- 
range attraction, derived from the zero-point fluctuations 
of the plasmon modes associated with two 2-D Wigner 
crystals |J. 

Recent attention has focused on fluctuation-induced 
forces in physical systems that contain charges on sur- 
faces. For example, correlation effects in some 2-D elec- 
tronic systems in semiconductor heterostructures |Io| l 
- specifically bilayer systems - could be viewed as a 
problem of charge-fluctuations on surfaces. Charge- 
fluctuation-induced forces at high temperatures may 
have another interesting realization in a collection of 
charged macroions in an aqueous solution of neutraliz- 
ing counterions, with or without added salt |Tl],[l^]. The 
macroions may be charged membranes, stiff polyelec- 
trolytes such as DNA, or charged spherical colloidal par- 
ticles. If the surface of macroions is highly charged, most 
of the counterions in solution condense onto the surfaces 
|p3| and their fluctuations may become important Jl4[ . 

Specifically, consider a planar surface with charge den- 



sity en, where e is the electronic charge. A neutralizing 
counterion in the solution experiences an (unscreened) 
electrostatic attractive force of magnitude ATrnlBksT, 
where Ib = ^ T ~ 7 A is the Bjerrum length for an 
aqueous solution of dielectric constant e = 80 (H2O), 
ks is the Boltzmann constant, and T is the tempera- 
ture. The length scale - the Gouy-Chapman length - 
at which the thermal energy balances the electrostatic 
energy, given by A = l/(47rig7i), defines a layer within 
which most of the counterions are confined. For suffi- 
ciently high charge densities A <C L, where L is the linear 
size of the macroion, the "condensed" counterions can be 
considered as a quasi-two-dimensional ideal gas of density 
n. To capture correlation effects at high temperature, it 
is sufficient to consider in-plane fluctuations about a uni- 
form charge distribution. Explicit calculations show that 
charge fluctuations lead to an attractive force between 
two such plates, which scales as d" 3 for large separations 
d |^2]| . This picture may explain the attractive interac- 
tion between two highly charged macroions, observed in 
experiments [[l5| and in simulations fl6|| . Note that this 
charge-fluctuation-induced force is similar in spirit to the 
van der Waals interaction, and indeed formally identical 
to the Casimir force between two partially transmitting 
mirrors at high temperatures fjjj ]. Although the avail- 
able experimental results are for high temperature, it is of 
fundamental interest to understand the low temperature 
interactions as well. Moreover, although it is unlikely 
to be relevant for macroions, such considerations may 
well impact on solid state systems, such as semiconduc- 
tor bilayers. The purpose of this Letter is to study this 
charge-fluctuation-induced force at zero-temperature. 

To this end, we consider a model system composed of 
two uniformly charged planes, separated by a distance d, 
each having a charge density of en. Confined to move 
on them are negative mobile charges, e.g. classical elec- 
trons, with particle density n on each of the surfaces, so 
that the system as a whole is neutral. At sufficiently low 
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temperature, the charges on the surface crystallize into a 
2-D triangular lattice as in Fig. [I]. This occurs when A is 
small compared to the average spacing between charges: 
X^/n <C 1, or r = e 2 ^/tui/ 'ekgT - the ratio of aver- 
age Coulomb energy among charges to thermal energy - 
is sufficiently large [ pj[ . For electrons on surface of liq- 
uid helium, it has been experimentally determined that 

r - 100 ra. 





FIG. 1. Two staggered Wigner crystals formed by the "con- 
densed" counterions. 

To estimate the attractive force arising from correla- 
tion effects at low T, we follow Ref. [^0| and consider 
first the limit of T — > oo, i.e., the ground state of the sys- 
tem. In this limit, it is easy to see that when two classical 
Wigner crystals experience each other's electric fields, the 
two lattices of charges stagger in order to minimize the 
energy of the system. In this case, the "classical" elec- 
trostatic force per unit area of such a charge pattern can 
be calculated, 
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where repeated indices are summed. Here, K a p(R — 

R) = <5r,R' Er» [ 0a/?(R " R") + A Q/3 (R " R" - c)] - 

<t> af 3(R - R), <j> a p(r) = d a dp-jff;r 7^ 0, and A Q/3 (r) = 

c> a 6^,3 J^ 2 =p . For a general lattice, the square of the 
phonon frequencies are the eigenvalues of the dynamical 
matrix |2lf| , which in this case may be written as 
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where m is the mass of the charges and K(k) and 
A(k) are 2x2 matrices whose elements are defined by 
K a/3 (k) = J2r K a f3(R) e" lk R and similarly for A Q/3 (k). 
Using an expansion in reciprocal lattice space, they are 
explicitly given by: 



K af3 (k) = A Q/3 (0) 
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for large d; where c is the relative displacement vec- 
tor between two lattices of the different plane. Hence 
there is a short-range attractive force, which decays ex- 
ponentially with the lattice constant as the characteristic 
length scale. In addition to this short-ranged force, we 
show below that zero-point fluctuations of the plasmon 
modes (charge fluctuations) also lead to an attractive 
long-range interaction, analogous to, but scaling differ- 
ently (in distance) from the standard zero-temperature 
Casimir effect. Moreover, at sufficiently low T, this force 
would still be manifested in the system of classical elec- 
trons. 

To study the zero-temperature charge-fluctuation- 
induced force, we first evaluate the phonon spectrum of 
the system of two coupled Wigner crystals. The positions 
of the charges are 
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where u A ' B (R) is the small deviation from the equilib- 
rium lattice sites R. Within the harmonic approxima- 
tion, the potential may be written as 
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where G are the reciprocal lattice vectors. In general, to 
obtain Wj(k), the frequency of the jth mode (j = 1, 4), 
we have to diagonalize D(k) numerically. This calcula- 
tion has been performed in Ref. [^2). However, since we 
are interested in the long-wavelength limit and large dis- 
tance asymptotics, we approximate D(k) in the following 
fashion: 
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where we have defined the matrix D°(k) with elements 



^(k) 



. This approximation entails neglecting 
contributions from higher order terms in k and from 
nonzero reciprocal lattice vectors, which are exponen- 



tially small (~ e ). Therefore, Eq 



is a good 
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approximation provided that d is larger than the average 
spacing between charges on the surface. 

Within this approximation, D(k) can be diagonalized 
to yield the following dispersion relations: 



jf(k) = 2A; uj 2 (k) = 2A 
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where we have chosen A(0) appropriate for staggered tri- 
angular lattices: A n = A 2 2 = A and A 12 = A 2 i = 0, 
with A ~ e~ Gd These modes can also be derived 

by treating the coupling between two isolated Wigner 
crystals as a perturbation. Mode 1 is one of the two 
optical modes which correspond to out-of-phase vibra- 
tions of the charges in opposite planes. The finite gap at 
k = vanishes exponentially with d for large distances as 
also found in Ref . |^] . Mode 3 is the transverse phonon 
mode, which describes the shear mode of the system, 
similar to that of a single 2D Wigner crystal. We remark 
that this approximation gives zero for its frequency, but 
upon including terms that involve next order in k the dis- 
persion is linear: u>3(k) ~ v s k p2| |. Its sound velocity v s 
is roughly a constant - the transverse sound velocity for 
an isolated Wigner crystal - plus a small correction (ex- 
ponentially decaying for large d) which arises from the 
interlayer coupling. In fact, all higher order terms are 
exponentially damped for large d. Modes 2 and 4 may 
be interpreted as the out-of-phase and in-phase plasmon 
modes, respectively. An interesting feature is that the 
out-of-phase plasmon mode has a gap at k = in the 
presence of the coupling. The in-phase plasmon mode 
vanishes as \fk as k — > and its sound velocity diverges. 
Physically, the transverse phonon and the in-phase plas- 
mon mode describe the charges in different planes os- 
cillating in-phase. We note that their existence has been 
shown to be a general property for a 2D Coulomb plasma 
p3| , and does not specifically depend on the nature of the 
underlying Wigner lattice. Thus, the results of this pa- 
per are universal, independent of the nature of the ground 
state. 

For the zero point energy, the dominant modes in 
the long wavelength limit are the plasmons: modes 2 
and 4. Neglecting exponentially small contributions, i.e. 
A — > 0, we obtain the zero-point energy (relative to the 
infinite separation) associated with the interactions be- 
tween the surfaces 
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where A is the area of the planar surface. This leads to 
an attractive pressure 
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where a is a positive numerical constant of order unity, 
explicitly given by 
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Thus, zero-point fluctuations induce a long-range at- 
traction which decays with a novel power law ~ d~ 7 1 2 . 
This should be contrasted with the usual Casimir-like 
force ~ d~ , which arises from, for example, the acoustic 
phonon zero-point fluctuations. We note that this power 
law stems from the 2 dimensional nature of charged sys- 
tems: 2-D plasmons do not have a finite gap, as they do 
in 3D. Note also that all the higher order terms in k, as 
well as those that we have neglected in our derivation, 
decay exponentially with d; therefore Eq. (|Io| ) is the 
dominant term for large distances. For an order of mag- 
nitude estimate, assuming m ~ 10 -25 kg, n ~ 1/50 A -2 , 
d ~ 10 A, and e - 80, we find n - 10~ 25 J/ A 3 . This is 
close to the magnitude of the "classical" force in Eq. (jl]) : 
F s /A ~ 10~ 24 J I A 3 , and thus may be just as important 
under suitable conditions. 

To recapitulate, we have argued that there is a long- 
range force at T — 0, derived from the zero-point fluc- 
tuations, which must be added to the zero-temperature 
"classical" force. At finite temperatures, an explicit cal- 
culation |25) using the Bose-Einstein distribution shows 
that at large separations, an additional contribution from 
the plasmon modes to the attractive force is of the form 
a2ksTd~ 3 , which agrees exactly (even the prefactor a.%) 
with the high temperatures result of Ref. |l2j . Moreover, 
the effect of finite T on the exponential force is to modify 
it with a "Debye- Waller" factor, weakening it, and even- 
tually causing it to vanish (25). Thus, for finite T, we 
have the following expression for the correlated attrac- 
tions for a system of two coupled planar Wigner crystals: 
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where ao,i,2 are constants and a is the range of the short- 
ranged attraction, of the order of the lattice constant. 

Finally, we comment that experimental observations 
of our result may prove subtle, as indicated by other ex- 
amples of phonon-fluctuation-induced interactions. For 
example, in Ref. [^6| the effect of zero-point fluctuations 
of phonons on the wetting transition of a He thin film 
was investigated, and the effect was shown to be small. 
Perhaps, a more appropriate system in which this be- 
havior may be manifested is bilayer quantum well sys- 
tems. Indeed, recent experimental techniques allow for 
2-D confinement of electrons, and a 2D plasmon disper- 
sion has been confirmed experimentally in GaAs/AlGaAs 
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heterostructures |27[ , and may exhibit fractional-power- 
law predicted in this Letter. Although there has been ex- 
perimental progress in observing single Wigner crystals 
PH , the study of interacting Wigner crystals remains an 
experimental challenge. In any event, the attractive force 
discussed in this Letter may assist with the conceptual 
understanding of correlation effects for highly charged 
surfaces and more fundamentally, indicates that for two 
neutral planes with mobile charges, the zero-point inter- 
action has a universal scaling different from that of the 
standard Casimir effect. 
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